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ABSTRACT 

We present results on the Soft Gamma Repeater (SGR) 0501+4516, discovered by the Swift Burst 
Alert Telescope (BAT) on 2008 August 22. More than 50 bursts were identified from this source, out 
of which 18 bursts had enough counts to carry out spectral analysis. We performed time-averaged 
spectral analysis on these 18 bursts using 8 models, among which the cut-off powerlaw and the two- 
blackbody models provided the best fit in the 15-150 keV energy range. The cut-off powerlaw model fit 
yields a mean photon index Tcpl = 0.54±0.11 and a cut-off energy Ec — 19.1±1.8 keV for the bursts. 
The mean hard and soft blackbody temperatures are found to be kTsB h — 12.8±0.7 keV and kTsB s 
= 4.6±0.5 keV, respectively, and are anti-correlated with the square of the radii of the hard and soft 
emitting regions {Rbb h an d Rbb s ) as R BB oc kT~ 5 8 and R 2 B b, kT~ 2 - 7 , respectively. The soft and 
hard component temperatures with different indices support the idea of two distinct emitting regions 
with the hard component corresponding to a smaller radius and the soft component corresponding to a 
larger radius, which further corroborate the idea of the propagation of extraordinary (E) and ordinary 
(O) mode photons across the photosphere, as predicted in the magnetar model. We notice strong 
burst fluence-duration correlation as well as hardness ratio-duration and hardness ratio-fiuence anti- 
correlations for the SGR 0501+4516 bursts. The burst fluences range from 



2.7x10 6 ergs cm 2 , consistent with those observed for typical short SGR bursts. 



4.4x10 9 ergs cm 2 to 
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1. INTRODUCTION 

The Swift mission (Gehrels et al. 2004), launched 
in November 2004 to explore the gamma-ray bursts 
(GRBs), has provided an excellent opportunity to de- 
tect and study the 7-ray activities from known Soft 
Gamma Repeaters (SGRs) as well as to discover new 
SGRs. SGRs, belonging to the class of 'magnetars' (Dun- 
can & Thompson 1992), are highly magnetized (B ~ 
1O14-15 G) and slowly rotating (P ~ 2—12 s) neutron 
stars characterized by short, bright bursts of hard X- 
rays and soft 7-rays (see Mereghetti 2008 for a recent 
review). During its active state, an SGR can go through 
periods of intense bursting activity lasting from a few 
days to months, however, it can also remain dormant for 
many years. The bursts, often varying in duration, are 
classified into three main categories: short, intermediate 
and giant flares. The short bursts are the typical kind 
observed during an outburst, marked by timescales ~ 
0.1-0.5 s and luminosities ~ 10 38 -10 41 ergs s _1 , whereas 
the intermediate bursts have timescales ~ 1-60 s and lu- 
minosities ~ 10 41 -10 43 ergs s _1 (Mereghetti 2008). Gi- 
ant flares, the rare and unique events, are distinguished 
by their extreme energies (10 44 -10 47 ergs s _1 ), long du- 
rations (~ hundreds of seconds), and the presence of a 
coherent pulsating decaying tail, consistent with the spin 
period of the neutron star. Persistent X-ray emission, in- 
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terpreted as originating due to the magnetospheric cur- 
rents driven by twists in the evolving ultra-high magnetic 
field (Thompson & Duncan 1995), has also been observed 
from SGRs in the 0.1-10 keV band with a typical X-ray 
luminosity of ~ 10 35 ergs s _1 , and the spectrum is gen- 
erally described by an absorbed powerlaw (PL) plus a 
blackbody component (Mereghetti 2008). 

According to the magnetar model (Thompson & Dun- 
can 1995), the dominant form of energy powering an SGR 
is its decaying ultra-strong magnetic field. The surface 
of the neutron star is heated and fractured by instabil- 
ities generating Alfven waves that accelerate electrons, 
and in turn give away their energy in short bursts. The 
model also suggests that the high-energy dissipated re- 
mains trapped in the magnetosphere as the 'trapped fire- 
ball', and it shrinks in size with time. Alternatively, the 
bursts can arise from heating of the corona by magnetic 
reconnection in the stellar magnetosphere resulting in 
intermediate- type flares (Lyubarsky 2002). The giant 
flares are caused by the sudden rearrangement of the 
star's magnetic field producing global crustal fractures 
(Thompson & Duncan 1995). 

SGR 0501+4516 is a recently discovered SGR by the 
Swift 7-ray observatory (Holland et al. 2008; Barthelmy 
et al. 2008). Rossi X-ray Timing Explorer (RXTE) ob- 
servations discovered a 5.7 s X-ray pulsating counterpart 
(Gogus et al. 2008), with a dipole magnetic field B = 
3xl0 14 G estimated from its period and spin-down rate 
(Woods et al. 2008), thus confirming the magnetar na- 
ture of the source. Soon after, the GLAST Burst Moni- 
tor (GBM) onboard Fermi satellite triggered and located 
the bursts from SGR 0501+4516 (van der Horst & Con- 
naughton 2008; McBreen et al. 2008). Multiwavelength 
observations reported the detection of its infrared (Tan- 
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TABLE 1 

Summary of SGR 0501+4516 bursts observed with Swift-BAT. 
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Note. — T100 is the total burst duration obtained by noting the start and the end time of the burst, manually. See j2]for details. 
a Total BAT exposure for the trigger numbers. 
b Count-rates are in the 15-150 keV energy range. 



vir et al. 2008; Rea et al. 2008) and optical (Fatkhullin 
et al. 2008; Rol et al. 2008) counterparts. However, no 
radio emission was detected from this source during the 
outburst (Kulkarni & Frail 2008; Gelfand et al. 2008). 
Suzaku observations of the SGR performed from 2008 
August 26-27 detected 32 bursts with the X-ray imaging 
spectrometer (XIS) and the hard X-ray detector (HXD). 
The persistent X-ray emission obtained with the XIS in 
the 0.4-10 keV range was best fitted by the combination 
of a blackbody (kT = 0.69±0.01 keV) plus a PL (F = 
2.8±0.1) component, modified by an interstellar absorp- 
tion N H = (0.89±0.08)xl0 22 cm- 2 (Enoto et al. 2009). 
The 1-200 keV combined XIS+HXD burst spectrum was 
described by a two-component blackbody model having 
soft and hard temperatures of kT — 3.31q'| keV and 
lS.lti'.g keV, respectively, and was interpreted as the 
population of ordinary (O) and extraordinary (E) mode 
photons propagating across the magnetosphere (Enoto 
et al. 2009). Konus-Wind 7-ray burst spectrometer ob- 
served the SGR bursts from 2008 August 23-26 and the 
20-200 keV spectra were best fitted by either a pow- 
erlaw with an exponential cut-off (CPL) model or an 
optically-thin thermal bremsstrahlung model (OTTB), 
both of which gave a peak energy (E p ) or tempera- 
ture (kT) in the range of 20-45 keV (Aptekar et al. 
2009). Subsequently, five XMM-Newton and two INTE- 
GRAL observations of the source were reported, with 
a hard X-ray variable source detected by INTEGRAL 
only during the first pointing (Rea et al. 2009). The 
phase-coherent timing analysis performed with XMM- 
Newton, Suzaku-XIS and Swift-XBT data further refined 
the period P = 5.7620695(1) s, period derivative P = 
6.7(1) xl0~ 12 s s~\ and magnetic field B ~ 7xl0 13 - 
2x 10 14 G, with the evidence of a second period derivative 
P = -1.6(4) xl0~ 19 s s~ 2 (Rea et al. 2009). Moreover, 
the spectrum of the outburst indicates a trend of spec- 
tral softening, with the blackbody component decaying 



slower than the PL component and the spectral evolution 
changing from phase to phase (Rea et al. 2009). 

In this paper, we report the Swift-BAT observations of 
the recently discovered magnctar SGR 0501+4516. The 
paper is organized as follows: ^and 321 describe the ob- 
servations and the burst spectroscopic analysis, respec- 
tively. In [JH we present the results of the time-averaged 
spectral and statistical analyses of the SGR bursts, which 
are then discussed in <JS] The conclusions are summarized 
in g3 

2. OBSERVATIONS AND DATA REDUCTION 

The initial outburst from SGR 0501+4516 was discov- 
ered with the Burst Alert Telescope (BAT) on 2008 Au- 
gust 22 (BAT calculated R.A = 05h 01m 04s, Dec = 
+45d 16' 20", with an uncertainty of 3' radius; Holland 
ct al. 2008; Barthelmy et al. 2008). Swift-BAT is a 
highly sensitive and large field of view (1.4 sr half-coded) 
hard X-ray telescope using a coded aperture mask oper- 
ating in the 15-350 keV energy range (Barthelmy et al. 
2005). The BAT detector plane is made of 32,768 pieces 
of CdZnTe (CZT: 4x4x2 mm) and the coded-aperture 
mask consists of ~ 52,000 lead tiles (5x5x1 mm) hav- 
ing a 1-m separation between mask and detector plane 
(Barthelmy et al. 2005). We used data from all the BAT 
triggers (see Table 1) and observation sequences (from 
00321174001 to 00321174061) available for this source. 
Intense bursting activity was observed between 2008 Au- 
gust 22-23 with the last BAT trigger on 2008 Septem- 
ber 3. The standard BAT software distributed within 
FTOOLS under the HEASoft packag<0 (version 6.4.1) 
and the latest calibration files available were used to an- 
alyze the data. The burst pipeline script, batgrbprod- 
mc0, was run to process the BAT trigger events. In this 
study, we screened out the faint bursts (characterized by 

^ http: / /heasarc. gsfc.nasa.gov/lheasoft/ 

® http://hoasarc.11asa.gov/docs/swift/a11alysis/throads/batgrbproductthrcad.html 
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a count-rate < 0.1 counts s _1 ) from the ~ 50 total bursts 
observed, leaving 18 bursts for the temporal and spectral 
analysis. In the cases where more than one burst was ob- 
served in a trigger, the bursts were denoted by the letters 
'A', 'B', etc., preceded by the trigger number (for exam- 
ple, the two bursts considered for spectral analysis from 
the trigger 000321174000 were named as 32 11 74 A and 
321174B). In Table 1, we present the summary of the 
bursts identified for spectral and temporal analysis along 
with their total BAT exposure and count-rates (counts 
s _1 ) in the 15-150 keV energy range. 

Background-subtracted light curves were made for all 
the bursts using the task batbinevt in the 15-25, 25-50, 
50-100 and 100-150 keV energy ranges binned at 4 ms 
time resolution after applying the maskweighting tech- 
nique using the batmaskwtevt task. Maskweighting is 
defined as the background-subtracted counts per fully il- 
luminated detector for an equivalent on-axis source and 
involves assigning each event a weight according to the 
illumination fraction of the detector it was detected in. 
A sample of the burst light curves is displayed in Figure 
1. To study the burst durations, we determined Tgo from 
the light curves by running the task battblocks, where Tgo, 
the standard parameter describing the burst duration of 
GRBs, is defined as the time to accumulate 5% to 95% 
of the observed photons. However, battblocks failed to 
determine the burst's duration for a few weak bursts. 
Hence, Xioo interval was used in our analysis, which was 
determined by noting the start and the end time of the 
burst emission, manually, and the Tioo values range from 
~ 0.03 s to ~ 0.5 s with an average value of 0.16±0.02 s. 
In Figure 2, we show the Tioo duration histogram of SGR 
0501+4516's bursts, overplotted with the Tgo duration 
histogram ranging from 0.07 s to 0.25 s, for comparison. 
The solid and the dotted curves are the lognormal Gaus- 
sian best-fits, peaking at 0.06±0.02 s and 0.09±0.02 s 
with a Gaussian width a of 0.8±0.1 s and 0.5±0.1 s for 
Tioo and Tgo, respectively. 

Response matrices for BAT spectra were generated for 
each of the 18 bursts using the batdrmgen task and the 
spectral fitting was restricted between 15-150 keV band, 
since the BAT mask becomes transparent around 150 
keV. Finally, a systematic error was applied to the BAT 
spectra using the batphasyserr task in order to account 
for residuals in the response matrix. The spectral analy- 
sis was performed using XSPEC vl2.4.1. Errors quoted 
are at the 90% confidence level. 

3. BURST SPECTROSCOPY 

The time-averaged spectral analysis was performed us- 
ing the Tioo interval. The bursts spectra were subse- 
quently fitted with 4 single component models: PL, CPL 
(characterized by an exponential cut-off energy Ec and a 
powerlaw photon index Tcpl), thermal bremsstrahlung 
(Brcmss), blackbody radiation with normalization pro- 
portional to the surface area (Bbodyrad in XSPEC; ab- 
breviated as BB throughout this paper), and 4 dou- 
ble models by adding a Bbodyrad component to all 
the above mentioned single models: BB+PL, BB+CPL, 
BB+Bremss and BB+BB. We first averaged the mea- 
sured total x 2 f° r a ll the bursts to determine their aver- 
age reduced \ 2 ((xt)-> where v is the number of degrees 
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TABLE 2 

Summary of the average \ 2 and F-test values. 
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Note. — Obtained for the 18 time-averaged burst spectra of 
SGR 0501+4516 using Swift- BAT. See fj3]for details. 



of freedom) and their F-test probability. The results, 
considered to give a qualitative analysis of the spectral 
model fit goodness, are summarized in Table 2 in the 
order of increasing (xl)- Among the single-component 
models, the CPL model gave the best fit ((xl) = 1.004) 
and a single PL model gave the worst-fit ((xt) — 2.766). 
The Bremss model, although considered best to describe 
the hard X-ray spectra of SGRs, also gave a large (xl) 
and failed to reproduce the spectral energy distribution 
for most of the bursts in the high-energy band, as re- 
ported by various authors in their study of other SGRs 
(for example, Feroci et al. 2004; Israel et al. 2008). 

We next investigated the significance of the double 
models: BB+PL, BB+CPL, BB+Bremss and BB+BB. 
As shown in Table 2, the BB+BB and BB+CPL models 
provide comparable (xl) values, but the BB+CPL model 
was ruled out because of the large F-test probability. For 
the BB+Bremss and BB+PL models, even though the F- 
test probability suggests significant improvement of the 
fit by the addition of a BB component, the (xt) values 
obtained are relatively high. Therefore, out of the 8 trial 
models, the CPL and BB+BB models yielded statisti- 
cally acceptable xl values. Figures 3 (a) and (b) show 
a sample of the best-fit spectra obtained with the CPL 
and BB+BB models for a burst (321583E) of duration 
- 0.33 s. 

4. RESULTS 

In the following subsections, we present the results 
from our spectroscopic analysis and also investigate the 
correlation between the spectral and temporal properties 
of the observed bursts from SGR 0501+4516. 

4.1. Time-averaged spectral parameters 

The best-fit time-averaged spectral parameters of the 
bursts are summarized in Table 3. We have deter- 
mined the cooler and hotter BB radii (Rbb 3 and Rbb h ) 
from the two-component BB model normalization K ~ 
R\ m /D\ , where Rk m is the radius of the emitting re- 
gion in km and Dio is the distance to the source in 
units of 10 kpc. Although the exact distance to SGR 
0501+4516 is not yet known, a distance of 1.5 kpc based 
on the proximity of its direction to the supernova rem- 
nant HB9 has been suggested by Gaensler & Chatterjee 
(2008). More recent studies used distance estimates of 
1.5 kpc (Aptekar et al. 2009), 5 kpc (Rea et al. 2009) 
and 10 kpc (Enoto et al. 2009). In this work, we as- 
sume a distance to the source in units of 10 kpc (£*io) 
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to be consistent with the XSPEC BB model normaliza- 
tion and scale all the derived parameters in units of D w . 
We determined a mean radius of ~ 7.2±1.5 Diq km and 
~ 0.9±0.2 Diq km for the soft and hard BB compo- 
nents, respectively, which are consistent with the cooler 

(~ S- 9 ^:? D w km ) and hotter (~ °- 46 ^o:i4 D w km ) 
BB radii obtained for SGR 0501+4516 using Suzaku ob- 
servations (Enoto et al. 2009). The histograms of the 
best-fit spectral parameters (Tcpl, Eq, kTsB 3 , kTBB h ) 
and the BB radii of the soft and hard components (Rbb 3 
and Rbb h ) are shown in Figure 4. In Table 4, we sum- 
marize the mean value obtained from a Gaussian fit to 
the histogram plots of all spectral parameters and the 
inferred BB radii (see Figure 4). The results obtained 
for SGR 0501+4516 are in good agreement with those 
obtained for the typical short-duration bursts observed 
in other SGRs such as SGR 1900+14 and SGR 1806-20 
(Olive et al. 2004; Feroci et al. 2004; Nakagawa et al. 
2007; Israel et al. 2008). 

The burst fluence, tabulated in Table 3, was esti- 
mated by multiplying the time-averaged flux (in units 
of ergs cm -2 s^ 1 obtained from the spectral fits) by its 
respective burst duration Tioo- For SGR 0501+4516, 
the fluence ranges from ~ 4.4xl0 -9 ergs cm~ 2 to 
~ 2.7xl0 -6 ergs cm -2 . Observational results sug- 
gest that SGR bursts fluence ranges from ^10~ 10 
~ 10~ 4 ergs cm -2 and follows a powerlaw distribution 
(Mereghetti 2008). The corresponding burst energies are 
in the range of ~ 4.2xl0 36 ergs to ~ 2.8xl0 39 ergs. As- 
suming isotropic emission, we calculate the burst peak 
luminosity as L = AttD^qF, where F is the time-averaged 
flux of the bursts (see Table 3), and find that it ranges 
from - l^j^xlO 38 D\ ergs s^ 1 to - 6.71^ xlO 40 
D\ Q ergs s _1 . We have also determined the luminosities 
corresponding to the soft and hard components of the 
BB+BB model; these are reported and discussed in [JSJ 

4.2. Statistical analysis of SGR 0501+4516 

Statistical studies have unveiled some basic properties 
such as the burst energy injection and radiation mech- 
anisms of the SGRs (Gogus et al. 1999, 2000, 2001). 
In this subsection, we consider the statistical properties 
of the SGR 0501+4516 bursts. First, we investigated 
the fluence distribution of the SGR bursts with duration 
and estimated the significance of this correlation using 
the Spearman test (Spearman 1904). The fluences show 
a strong positive correlation with Tioo, with a PL index 
of 2.0+0.2 as shown in Figure 5. The Spearman rank 
order test applied to SGR 0501+4516 burst fluences and 
durations yielded a correlation coefficient p = 0.9 and the 
probability (P) that this correlation is due to a random 
fluctuation is 3.0xl0 -7 corresponding to 95% confidence 
level. Such a correlation has been observed by Gogus et 
al. (2001) for SGR 1806-20 and SGR 1900+14, which 
was interpreted to be similar to the PL relation between 
the total energy and duration of the earthquakes estab- 
lished by Gutenberg & Richter (1956). 

Next, we examined the burst spectral variations with 
temporal properties to draw comparison to other SGRs. 
For that, we computed the hardness ratio, defined as 
HR = (H-S)/(H + S), where H and S are the back- 
ground subtracted hard and soft photon counts in the 
25-150 keV and 15-25 keV energy bands, respectively. 



Figures 6 (a) and (b) show HR plotted against the burst 
duration and fluence, respectively. We clearly see an 
anti-correlation in both cases and a powerlaw fitted to 
the data gave indices of -0.20+0.03 and -0.10+0.02, re- 
spectively. This anti-correlation was further quantified 
by carrying out the Spearman test and we obtain a cor- 
relation coefficient of p = -0.7 (P = 6.5xl0~ 4 ) for HR 
vs. burst durations and p = —0.8 (P = 1.8xl0 -4 ) for 
HR vs. burst fluences. We discuss the observed correla- 
tions in ij5] a detailed analysis on the spectral evolution 
in SGR 0501+4516 will be presented elsewhere. 

5. DISCUSSION 

We have carried out a comprehensive study of the spec- 
tral and temporal properties of the recently discovered 
SGR 0501+4516. The observed burst durations in SGR 
0501+4516 typically range from ~ 0.03-0.5 s and fol- 
low a lognormal distribution as seen in the case of other 
SGRs (for example, Gogus et al. 2001; Woods et al. 
1999). T 9 o durations of ~ 0.093 s and - 0.162 s were 
measured for SGRs 1806-20 and 1900+14 (Gogus et al. 
2001), respectively, and ~ 0.099 s for the anomalous X- 
ray pulsar (AXP) IE 2259+486 (Gavriil et al. 2004). 
The Tioo duration of the bursts observed from SGR 
0501+4516 showed a lognormal distribution peaking at 
~ 0.06+0.02 s with an average value of ~ 0.16+0.02 s 
for the 18 bursts. The short bursts observed from SGR 
0501+4516 using Swift-BAT were fitted with 8 spectral 
models, out of which the single-component CPL and the 
two-component BB+BB models provided the best fit. 
Feroci et al. (1994) also suggested CPL and BB+BB as 
best-fit models for 10 short bursts from SGR 1900+14 ob- 
served using BeppoSAX in the 1.5-100 keV energy range. 
We have also explored the possible correlations between 
the spectral and the temporal burst properties in SGR 
0501+4516 and the implications of these observed fea- 
tures are discussed below in the context of the magnetar 
model (Thompson & Duncan 1995). 

As shown in gJJ the SGR 0501+4516 burst fluence 
follows a powerlaw distribution (a = 2.0+0.2; see Fig- 
ure 5) with the burst duration. Similar correlations have 
been found for earthquakes with PL indices ranging from 
1.4 to 1.8 (Gutenberg & Richter 1956; Chen et al. 1991) 
and for solar flares with a PL index ranging from 1.53 to 
1.73 (Crosby et al. 1993). Earthquake-like behaviour for 
SGRs were first pointed out by Cheng et al. (1996) using 
SGR 1806-20 with a PL index of 1.66. The PL relation 
between the seismic moment (oc energy) and duration 
was shown by Lay & Wallace (1995), which yielded an 
index of 3.03. Gogus et al. (2001) has also successfully 
investigated this correlation and obtained a PL index 
of 1.05+0.16 and 0.91+0.07 for SGR 1806-20 and SGR 
1900+14, respectively. This behaviour shown by earth- 
quakes was also interpreted on the basis of self-organized 
criticality (Bak et al. 1988), providing a framework for 
understanding complexity and scale invariance in sys- 
tems showing irregular fluctuations. Self-organized criti- 
cality can be considered as a characteristic state of criti- 
cality formed by self-organization in a long transient pe- 
riod at the border of stability and chaos with the events 
following a powerlaw distribution. However, this theory 
cannot predict the strength or time of the next event 
(Bak et al. 1988). Moreover, we notice that this is the 
first time the PL index obtained from an SGR matches 
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TABLE 3 

Summary of the spectral fit results to the SGR 0501+4516 bursts 
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Note. — The radius, fluence, flux and luminosity values quoted are in the 15-150 keV energy range, assuming a distance to the source 
in units of 10 kpc (Z?io). Errors quoted are at the 90 % confidence level. See [[4] for details. 



TABLE 4 

Summary of the Gaussian fit to the spectral parameters. 



Parameters 


Fcpl 


E C 


kT BBh 
(keV) 


kT BBs 


R BB h 


Rbb s 






(keV) 


(keV) 


(km) 


(km) 


Gaussian Mean of parameters 


0.54±0.11 


19.1±1.8 


12.8±0.7 


4.6±0.5 


0.9±0.2 


7.2±1.5 


Width of Gaussian distribution (tr) 


0.38±0.11 


6.7±1.8 


1.2±0.9 


1.7±0.6 


0.6±0.2 


5.9±1.5 



quite well with those measured in both earthquakes and 
solar flares, suggesting that SGR bursts are analogous to 
earthquakes and have both crustal and magnetosphcric 
origins. 

We have also investigated the hardness-fluence and 
hardness-duration correlations for the 18 bursts observed 
from SGR 0501+4516. As suggested by other authors, 
hardness-duration and hardness-fluence anti-correlations 
are observed, as shown in Figures 6 (a) and (b), respec- 
tively. Moreover, as pointed out by Gogus et al. (2001), 
we see that the SGR bursts soften with increasing burst 
durations, the details of which will be presented in a 
follow-up paper. The hardness-fluence anti-correlation 
(see Figure 6 (b)) has been explained by Gogus et al. 
(2001) as due to the emitting plasma in local thermody- 
namic equilibrium with the radiative area decreasing at 
lower fluences, and secondly, due to the spectral intensity 
of the radiation field being below that of a BB, causing 
the emitting plasma temperature to remain in a narrow 
range and higher at lower luminosities. However, it is 
also suggested that these two possibilities depend on the 
rate of energy injection into the atmosphere (Gogus et 
al. 2001). 

The cut-off energy Ec obtained from the CPL model 
peaks around — 19.1±1.8 keV (Figure 4(b) and Table 4), 
consistent with that obtained using Konus- Wind obser- 
vations {Ec in the range of 20-45 keV; Aptekar et al. 
2009). The Tcpl for SGR 0501+4516 is also in good 
agreement with those obtained for SGR 1900+14 and 
SGR 1806-20 (see Table 3 of Feroci et al. 2004 and 



Tables 9 & 10 of Nakagawa et al. 2007). The mag- 
netar model (Thompson & Duncan 1996) suggests that 
the nonthermal emission observed in the SGRs origi- 
nate from the hydromagnetic winds of particles in the 
magnetosphere and the cut-off energy Ec can be asso- 
ciated with the plasma energy distribution in the mag- 
netosphere. This contributes to the hard X-ray emission 
through particle acceleration leading to comptonization 
and particle bombardment of the surface (Thompson & 
Duncan 1996; Mereghetti 2008) and/or from pairs cre- 
ated higher (~100 km) in the magnetosphere (Thompson 
& Beloborodov 2005). 

The spectra of SGR 0501+4516 were also well-fitted 
with a two-component BB model, with the low- and 
high-BB temperature clustered around 4.6+0.5 keV and 
12.8+0.7 keV, respectively (Figure 4 and Table 4). These 
values are consistent with those reported by Enoto et al. 
(2009) using Suzaku observations (kTgB s = 3.3t u '4 keV 
and kT BBh = 15.1±^;| keV). The soft and hard BB com- 
ponent gave a Gaussian mean radii Rbb 3 = 7.2+1.5 Diq 
km and RbBh = 0.9+0.2 D w km, respectively. Sev- 
eral lines of interpretation were made to account for the 
BB+BB emission. Many recent studies have shown that 
the broadband (1-100 keV) spectroscopic studies of short 
SGR bursts are well approximated by the sum of two 
blackbodies. Olive et al. (2004) reported that the time- 
integrated 2-150 keV energy spectrum of an interme- 
diate flare from SGR 1900+14 observed with HETE-2 
was best described by the sum of two BBs with tern- 
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peratures ~ 4.3 keV and ~ 9.8 keV. These results were 
in line with those obtained for the short bursts from 
SGRs 1900+14 (kT ~ 3.23 keV and - 9.65 keV) and 
1806-20 (kT ~ 4 keV and - 11 keV) by Feroci et al. 
(2004) and Nakagawa et al. (2007). Similar results were 
also reported for the burst spectra of SGR 1900+14 and 
SGR 1627—41 with temperatures kT BBs ~ 3—5 keV and 
kT BBh ~ 9-10 keV (Israel et al. 2008; Esposito et al. 
2008). Olive et al. (2004) interpreted the high-BB com- 
ponent to be arising from a trapped fireball, with the 
high-BB temperature being consistent with the theoret- 
ically predicted temperature of a trapped fireball (^11 
keV) and the low-temperature BB with the constant ra- 
dius (almost 30 times larger than the average radius of 
the fireball) as originating from the surface of the neu- 
tron star. Furthermore, they suggest that the emission 
regions observed may be due to the radiation transfer ef- 
fects in a superstrong magnetic field. In addition, these 
two BB temperatures can be considered as arising from 
the propagation of ordinary (O) and extraordinary (E) 
photons across the magnetosphere, with the photosphere 
of the E-mode photons located closer to the neutron star 
surface and the scattering photosphere of the O-mode 
photons higher up in the magnetosphere as suggested by 
Israel et al. (2008). Our results are also consistent with 
these interpretations based on the magnetar model. 

In Figure 7, we have plotted the soft and hard com- 
ponent BB temperatures as a function of the square of 
their corresponding radii to investigate any possible cor- 
relations. The size of the emitting regions derived from 
the soft BB temperature, though does not appear to have 
a constant value, approximates the expected radius of 
the neutron star and a powerlaw fit to the data yields 
an index of —2.7+1.2. Similarly, we fitted a powerlaw 
to the hard BB temperature and their respective radii, 
and we obtain a PL index of —5.8+1.0. The Spearman 
correlation test gives the following anti-correlation re- 
sults: p = —0.5 (P = 0.06) for the soft BB temperature 
and p = —0.7 (P = 0.01) for the hard BB temperature. 
This clearly demonstrates the existence of a strong anti- 
correlation between kTh and R 2 BBh ■ Olive et al. (2004) 
found a nearly constant radius for the low-BB tempera- 
ture, independent of the temperature whereas the high- 
BB temperature showed a clear anti-correlation between 
the radius and the temperature (see Figure 7 of Olive 
et al. 2004). Such a correlation (R 2 oc fcT~ 3 ) has been 
observed for SGR 1900+14 by Israel et al. (2008) during 
an intermediate outburst for the luminous phases of the 
flares (L tot > 10 41 ergs s" 1 ). Israel et al. (2008) iden- 
tified the presence of a natural separation region around 
25-30 km which corresponds to a critical surface where B 
= Bqed, as predicted in the magnetar model (Duncan 
& Thompson 1995). However, from our data on SGR 
0501+4516, we notice that the minimum radius of the 
soft BB overlaps with the maximum radius of the hard 
BB component, thereby suggesting that a fraction of the 
E-mode photons can be locally converted to the O-mode 
by Compton scattering and by photon splitting if the ef- 
fective temperature is high enough (Duncan & Thomp- 
son 1995). 

We have determined the soft and hard BB luminosities 
(L BBs and L BBh ), and plotted them against each other 
to explore their correlation. Figure 8 shows the bolomet- 



ric luminosities of the two BBs for the bursts fitted with 
the BB+BB model. While Feroci et al. (1994) found 
a constant ratio for the bolometric luminosities of the 
bursts, we do not. However the quantities are well corre- 
lated with a correlation factor p = 0.7 (P = 6.1 xlO -3 ). 
We also fit the data with a powerlaw and obtain an in- 
dex of 0.6+0.1, which is consistent with those obtained 
by Israel et al. (2008) for SGR 1900+14. The lumi- 
nosities obtained for SGR 0501+4516 fall below - 10 41 
ergs s _1 , and hence, with the current statistics, we can- 
not make an inference on the saturation effect of the low 
BB luminosity as seen for SGR 1900+14 (Israel et al. 
2008). The maximum bolometric luminosity obtained for 
the SGR 0501+4516 bursts is for the soft BB component 
with L BBs = 4.1 x 10 40 D\ ergs s _1 corresponding to the 
radius R BBs = 14.2 D w km and temperature kT BBs = 
6.7 keV. We compare this with the magnetic Eddington 
luminosity L B dd,B (Paczynski 1992; Thompson & Dun- 
can 1995) given by 

L Edd . B ~ 2 x 10 40 (-?-) 4 / 3 (-^-) 2 / 3 ergs s" 1 (1) 
OQED Mns 

This gives a B-field value of - 6.3t^xl0 13 G, which is 
consistent with the magnetic field strength of the dipolar 
component in the range ~ 7x 10 13 < Bd< 2x 10 14 G (Rea 
et al. 2009). 

6. SUMMARY AND CONCLUSIONS 

We have presented a quantitative analysis of the newly 
discovered SGR 0501+4516 observed by Swift-BAT . The 
spectra of the bursts from SGR 0501+4516 in the 15- 
150 keV energy range are best-fitted by both the CPL 
and BB+BB models. In the BB+BB model, the hot- 
ter temperature represents the smaller trapped fireball 
regions and the colder temperature corresponds to the 
regions consistent more or less with the radius of the 
neutron star. These temperatures are also associated 
with the propagation of E-mode and O-mode photons 
across the photosphere as predicted in the magnetar 
model. The nonthermal emission observed from the short 
bursts seems to originate from the magnetosphere. All 
the above findings are consistent with the typical short 
bursts observed from other SGRs and with the magne- 
tar model prediction of SGRs. However, we believe that 
more burst observations will help in further constraining 
the SGR emission mechanisms. 
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Fig. 1. — Examples of the background-subtracted 4 ms lightcurves of SGR 0501+4516 bursts detected by Siui/f-BAT in the 4 energy 
channels (Black: 15-25 keV, Red: 25-50 keV, Blue: 50-100 keV, Green: 100-150 keV) for Trigger numbers 00321481000 (Burst: 321481A; 
left) and 00321551000 (Burst: 321551; right). 
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Fig. 2. — Distribution of the Tioo and Tgrj durations of the SGR 0501+4516 bursts. The solid curve is obtained by fitting a lognormal 
Gaussian model, which peaks at 0.06+0.02 s. Also shown by a dotted curve is the lognormal fit to the Tgo durations, peaking at 0.09+0.02 s. 
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Fig. 3.— Sample spectra of an SGR 0501+4516 burst (321583E) of duration Tioo = 0.33 s. (a) CPL model fit with {xl) = 1.250 (56). 
(b) BB+BB model fit with (xl) = 1.009 (55). 
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Fig. 4. — Histogram plots for the best-fit spectral parameters fitted with a Gaussian. The spectral indices (Tqpl), and cut-off energy 
Ec are obtained from the CPL model. The soft BB temperature kTsB s > hard BB temperature kTgg h , and the inferred soft and hard BB 
radii (Rbb s ; RflB h ) are obtained from the BB+BB model. Details of the Gaussian fit are given in Table 4. 
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Fig. 6.— (a) Hardness ratio vs. Tioo fitted with a PL of index -0.20±0.03 (p = -0.7, P = 6.5xl0" 4 ). (b) Hardness ratio vs. fluence 
fitted with a PL of index -0.10±0.02 (p = -0.8, P = 1.8X1Q- 4 ). 
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Fig. 7. — Square of the radii of the BB+BB model components as a function of their temperatures. The red circles and blue squares 
denote the soft and hard kT values, respectively. The soft kT values show an anti-correlation of p = —0.5, P = 0.06 with R BB , and the 

hard kT shows a strong anti-correlation of p = —0.7, P = 0.01 with R BB ■ 
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Fig. 8. — Bolometric luminosity of the hard BB component (L BBh ) vs. that of the soft BB component (L BBg ) in units of 10 39 ergs s 1 . 
The solid line is the powerlaw fit which gives the relation L BBs = (L BBh ) a where a = 0.6±0.1. 



